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ABSTRACT 

In this paper we show that perturbations of the accretion flow within the central engines of some 
active galactic nuclei (AGNs) are likely to form shock waves in the accreting plasma. Such shocks, which 
may be either collisional or collisionless, can contribute to the observed high-energy temporal and spec- 
tral variability. Our rationale is the following: Observations show that the continuum emission probably 
originates in an optically thin, hot plasma in the AGN central engine. The flux and spectrum from this 
hot plasma varies significantly over light crossing timescales. Several authors have suggested that macro- 
scopic perturbations contained within this plasma are the sources of this variability. In order to produce 
the observed emission the perturbations must be radiatively coupled with the optically thin hot matter 
and must also move with high velocities. We suggest that shocks, which can be very effective in 
randomizing the bulk motion of the perturbations, are responsible for this coupling. Shocks should 
form in the central engine, because the temperatures and magnetic fields are probably reduced below 
their virial values by radiative dissipation. Perturbations moving at Keplerian speeds, or strong non- 
linear excitations, result in supersonic and super-Alfvenic velocities leading to shock waves within the 
hot plasma. 

We show that even a perturbation smaller than the emitting region can form a shock that significantly 
modifies the continuum emission in an AGN, and that the spectral and temporal variability from such a 
shock generally resembles those of radio-quiet AGNs. As an example, the shock inducing perturbation in 
our model is a small main-sequence star, the capturing and eventual accretion of which are known to be 
a plausible process. We argue that shocks in the central engine may also provide a natural triggering 
mechanism for the “cold” component of Guilbert & Rees two-phase medium and an efficient mecha- 
nism for angular momentum transfer. Current and future missions, such as ASCA, XTE, XMM, AXAF, 
and ASTRO-E may determine the importance of shock-related emission from the central engines of 
AGNs. 

Subject headings: galaxies: nuclei — galaxies: Seyfert — quasars: general — shock waves 


1. INTRODUCTION 

Recent observations indicate that the continuum emis- 
sion from AGNs originates within a hot plasma (e.g., 
Mushotzky, Done, & Pounds 1993, hereafter MDP93, and 
references therein). Variability arguments suggest that the 
plasma is contained within a region of size R x ~ t d c, where 
t d is the flux doubling timescale. Within the commonly 
accepted accretion model for AGNs, the hot plasma is 
probably located in the vicinity of a central supermassive 
black hole (Rees 1984, and references therein; Blandford 
1990, hereafter B90). In this model, most of the continuum 
emission originates from within about twenty gravitational 
radii of the black hole in the region generally referred to as 
the AGN central engine. 

An extensively considered accretion model of AGNs, the 
simple optically thick geometrically thin disk model that 
does not include hot plasma, does not adequately explain 
the continuum UV-to-y ray spectrum from AGNs. Several 
modifications to the thin disk geometry have been pro- 
posed: (1) the geometrically thick accretion disk supported 
by radiation pressure (Abramowicz, Calvani, & Nobili 
1980) with a geometry that is susceptible to global pertur- 
bations (Papaloizou & Pringle 1984), (2) the geometrically 
thick, optically thin disk model with two-temperature rela- 
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tivistic plasmas, first suggested for an X-ray binary 
(Shapiro, Lightman, & Eardley 1976), but which has been 
extended to AGN counterparts (e.g., Rees et al. 1982; White 
& Lightman 1989, hereafter WL89; Tritz & Tsuruta 1989; 
Narayan & Yi 1994), and (3) a disk-corona model (Liang & 
Price 1977; Haardt & Maraschi 1991; Zycki, Collin- 
Souffrin, & Czerny 1995, hereafter ZCC95; Tsuruta & 
Kellen 1995, hereafter TK95). For the two-temperature 
torus model Narayan & Yi (1994) propose that at lower 
accretion rates advection lowers the efficiency of the accre- 
tion process and stabilizes the otherwise unstable hot elec- 
trons, which suggests that this model may be more 
appropriate for sub-Eddington sources (see also Artemova 
et al. 1996). A common feature of all models that attempt to 
reconstruct the higher energy (~ 1 keV 1 MeV) power-law 
spectrum is the introduction of an optically thin hot plasma 
(with the electron temperature T e > 10 8 K) that may domi- 
nate the emission from the central engine. 

Large-amplitude temporal variability of the hot plasma 
emission is present in at least 50% of radio-quiet AGNs, 
and significant spectral variability of similar timescales is 
also very common (MDP93, and references therein; Green 
1993; McHardy 1988). The temporal and spectral variabil- 
ity has been generally explained in the context of pertur- 
bations on the thin, cold disk model, but rarely in the 
context of the geometrical models mentioned above. 
Models for the variability of radio-quiet AGNs include, for 
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example, the hot spots model (Abramowicz 1990; Wiita 

1993, hereafter W93), where a multitude of randomly dis- 
tributed flares or “ hot spots ” embedded within a thin disk 
are responsible for the total variability. Other models also 
assume randomly distributed perturbations on the thin disk 
(Pudritz & Fahlman 1982; Devries & Kuijpers 1989). The 
temperature of the thin disk in the central engine is found to 
be of order T ~ 10 6 K giving a sound speed of order c s ~ 
10 ~ 2 c or less (Frank, King, & Raine 1992, hereafter 
FKR92). The total volume affected by a single perturbation 
is therefore too small to easily account for the doubling of 
the entire continuum luminosity, making these models 
dependent on the existence of many perturbations of char- 
acteristic size, time, and brightness occurring on the disk 
surface. In contrast, within the hot plasma region, a dis- 
turbance may propagate at a much greater speeds, affecting 
a larger volume of the continuum emitting material. A 
single perturbation effectively coupled to the hot plasma 
may transfer enough energy to the radiating region within 
typical flux doubling timescales to account for the observed 
continuum variability. 

Shocks have been investigated in the context of radio- 
loud AGNs with shocks propagating down a jet and inter- 
acting with irregularities in the jet material (Qian et al. 
1991), or other forms of shocks (W93). In this paper we 
investigate the possibility that a strong perturbation, 
moving at supersonic or super-Alfvenic speeds in the hot 
plasma of the AGN central engine, forms a shock that 
modifies the continuum emission (see also Sivron & 
Tsuruta 1994). Such shocks may effectively couple macro- 
scopic perturbations with the optically thin hot matter, 
modifying the continuum spectrum. Shocks are also effec- 
tive in randomizing the bulk motion of the perturbations, 
which may have important consequences for the outward 
transfer of angular momentum. 

In our work, we utilize recent studies of the effect of pairs 
on the central engine structure (Guilbert & Rees 1988, here- 
after GR88; Lightman & White 1988, hereafter LW88; 
Coppi & Blandford 1990, hereafter CB90; Sivron & 
Tsuruta 1993, hereafter ST93; and Ghisellini & Haardt 

1994, hereafter GH94). These authors suggest that in 
compact central engines, dampening processes readily cool 
a significant fraction of the hot matter, radiating the excess 
energy away. The cooled matter then forms a cold “ phase ” 
component. We show that in central engines with near- 
Eddington luminosity the postshock matter is, in effect, an 
“increased compactness” region in which the shocked 
matter is cooled. The shocks are thus a natural mechanism 
by which radiative power is increased and cold phase 
matter is created. The cooling of postshock matter results in 
radiation from the shock front, which may be responsible 
for the “ soft ” flares, whereas cold phase matter behind the 
shock front may be responsible for the relatively hard 
“ dips ” in the light curve (see § 4). 

For mathematical simplicity we use a model in which the 
source of the perturbations is local, for example, a small 
main-sequence star. Such a star can be captured and even- 
tually accreted into the central engine (Syer, Clarke, & Rees 
1990, hereafter SCR90). For a wide range of accretion 
parameters, such an event naturally results in the creation 
of shocks. Such shocks can significantly modify the light 
curve from the X-ray emitting region of AGNs. 

In § 2 we show that macroscopic inhomogeneities, which 
we assume travel at Keplerian speeds, naturally posess 


supersonic or super-Alfvenic speeds when present in the hot 
plasma of the AGN central engine for the three models 
discussed above. In § 3 we show that the shocks effectively 
couple the perturbations with the hot plasma, and large 
enough perturbations can therefore result in the modifi- 
cation of the total emission. In § 4 we calculate the effects of 
a large shock on the light curve of a model AGN. The main 
purpose of this exercise is to demonstrate the type of spec- 
tral variability expected when shocks modify the emission 
from the central continuum source. We discuss our results 
in § 5, and present our conclusions in § 6. 


2. SHOCK FORMATION 

Recent observations as well as theoretical analyses 
suggest that hot plasma must exist in the central engine. 
From an observational point of view, the variable contin- 
uum of nearly power-law shape in the X-ray to y-ray bands 
could be explained as emission from hot plasma radiating 
synchrotron self-Compton or Comptonized bremsstrahlung 
radiation from an accretion disk, filaments, or clouds (B90; 
ST93; GH94; TK95). Theory predicts that the simple thin 
accretion disk will be disrupted within a critical radius 
where heating due to gravitational and viscous energy 
release surpasses the cooling by bremsstrahlung, the tem- 
perature rises and the gas pressure disrupts the thin disk 
(Rees et al. 1982; B90). Other plasma instabilities may also 
disrupt the thin disk (see, e.g., W93). The accretion flow 
inside this radius is not a simple thin disk and it may form a 
geometrically thick torus (cases [1] and [2] above), or a 
disk corona structure with dominant radiation coming from 
corona (Haardt & Maraschi 1991). 

In order to describe the physical conditions within the 
central engine in the above models we use the following 
parameters: The dimensionless radial distance to the super- 
massive black hole is r = R/R where R g = GM/c 2 is the 
gravitational radius. The bulk of X-ray continuum is 
thought to originate from within R x ~ 20 R g . The central 
black hole mass, A/, for Seyfert galaxies and radio-quiet 
quasars has been estimated from X-ray time variability data 
to be M > ~ 10 6 A/ g , and we therefore write the black hole 
mass in terms of M 6 = M/10 6 M 0 . The accretion rate, M, 
has been related to the X-ray luminosity, L x , through M = 
Lx/fac 2 ), where y\ is the efficiency of converting gravitational 
energy of the accreted matter into radiation. Luminosities 
of order L x — io 4I -*4 er g s -i are typical for the 2-10 keV 
band. The dimensionless accretion rate is m = M/M E , 
where 


M e = 


4nGMm p 
(< t t c) 


( 1 ) 


is the Eddington accretion rate, m p is the proton mass, and 
<r x is the Thompson cross section. Assuming tj ~ 0.1 and 
M 6 ~ 1 we find typical accretion rates of order 0.01 < 
~m < ~ 10. 


2.1. Collisional Shocks 

That collisional shocks form and Rankine-Hugoniot rela- 
tions can be written, when wave steepening occurs in rela- 
tivistic thermal plasmas was shown by Taub (1948). 
Collision-dominated waves and shocks may be present 
when the particle deflection length is smaller than the size of 
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a typical large perturbation in the central engine, i.e., when 


2(kTf 
Tin, e 4 In A 


< Rg , 


(2) 


where n ( is the ion number density, T e is the electron tem- 
perature, and the Coulomb logarithm is here taken to be 
In A ~ 20 (see, e.g., FKR92). Assuming quasi-spherical acc- 
retion geometry with solid angle fl, we may easily relate the 
disk density to the accretion rate: 

M = QR 2 nm p n ( r acc . (3) 

Here nm p is the average ion mass and v acQ is the radial 
accretion velocity. Equations (2) and (3) then give a lower 
limit to the accretion rate for AGNs that may support col- 
lisional shocks: 


m > 0.08 



where & e = kTJm e c 2 is the dimensionless electron tem- 
perature and ©„ = l/r is the dimensionless virial tem- 
perature. One can see that collision-dominated waves may 
form in the near or super-Eddington regime. This result is 
hardly surprising, since it is related to that of GR88, which 
claims that the central engine is almost always optically 
thick for near Eddington luminosity, with the main differ- 
ence being that the mean free path for photons in GR88 is 
replaced by the deflection length for electron collisions. We 
can compare the above result with the condition in Rees et 
al. (1982), m cr < 50a 2 on the maximal accretion rate such 
that ions and electrons decouple to form a two-temperature 
ion torus. Here a is similar to the parameter in the standard 
a accretion model, which relates the stress and pressure in a 
thin disk. We see that even two-temperature tori can mar- 
ginally support collisional waves. 

For the cases in which the shocks are collision dominated 
we use the thermodynamic parameters, such as the adia- 
batic constant T, the speed of sound c s , and others, from a 
paper on shocks in relativistic collisional plasmas (Iwamoto 
1989, hereafter 189, and references therein). The speed of 
sound in the plasma is 



1 

i + c 2 {k5(o/k 2 (o - mmM 



X 


2-i 


El Q ]- 1 

k 2 ( 0 . 


(5) 


where £ = mc 2 /kT, and Kj© is the modified Bessel func- 
tion of the second kind, which naturally arises from the 
inclusion of relativistic energies in the exponent of the parti- 
tion function (189, eq. [A22]). 

The equilibrium temperature does not quite reach its 
virial value because the outward transfer of angular 
momentum, which also gives rise to radiation, is dissipative. 
We therefore write the dimensionless ion temperature, 
©. = kTJnm p c 2 , as a fraction S of the virial temperature, 

0 , = 5 & v , ( 6 ) 


where S ~ < 1. 

The ideal gas law, which holds for collisional plasmas at 
high temperatures (189, eq. [A8]), is a reasonable approx- 
imation for the equation of state, and in the ion-dominated 
case the speed of sound in equation (5) reduces to 

c s ~ v ' / T^r _1/2 c , (7) 


where we assumed 0, < 1. If the electrons in the preshock 
region are coupled with the ions, the electron speed of 
sound is calculated using the exact form of equation (5). 
This yields a speed that is roughly c 5e ~ c/3, because the 
electrons are relativistic and may be described by a simple 
equation of state as p e ~ € e /3, where p e and e e are the pres- 
sure and the energy, respectively. Cases in which the pre- 
shock electrons are decoupled from the ions are not 
described in 189 or other references, and are therefore not 
discussed in this paper. We note, however, that the limit in 
equation (10) at the end of this section is still valid. 

An object or perturbation in the accreting plasma that is 
more compact than the surrounding material, so that it 
cannot thermalize or lose its angular momentum on 
dynamical timescales, moves at a velocity comparable to its 
Keplerian or free-fall velocity (see also § 3.1 and the next 
paper in this series, Sivron, Tsuruta, & Caditz 1996b, here- 
after STC96). The relative speed of such a compact object is 
therefore 

0 * = |t>* - v g \ = br~ l{2 c , (8) 

where v * is the velocity of the object, v g is the velocity of the 
accreting matter and b is a geometrical factor of order unity. 
The Mach number is therefore given by 

^ = r 1 / 2 . (9) 

C, 

This Mach number is insensitive to the object’s distance 
from the black hole because both the speed of sound and 
the typical velocity are approximately proportional to the 
square root of the virial energy. For a wide range of accre- 
tion rates 0.003 < 5 < 1 (189; WL89), which results in 

0.8 < Jt < 16 , (10) 

where for the upper limit we substitute T = 4/3. This result 
suggests that shock waves are likely to form in AGN central 
engines that satisfy the constraint imposed by equation (4). 
The Mach surface opening angle, 4> ~ 2 arcsin (c s ,/r t ), may 
be significant (see § 4). The shocked surface formed behind a 
small perturbation roughly takes a conical shape that is 
likely modified by the presence of strong gravitational and 
shear forces, radiation pressure, turbulent flows, and density 
gradients within the central engine. 

2.2. Collisionless Shocks 

For accretion rates lower than the limit of equation (4) 
the deflection length is larger than the characteristic size of 
the central engine, and collision-dominated shocks do not 
form. However, magnetic and electric fields and turbulence 
are expected to be present and collisionless shocks may 
occur (Kennel, Edmiston, & Hada 1985, hereafter KEH85). 
In this case, the plasma state across the shock is determined 
not by the binary collisions, but by collective interactions 
between particles and self-consistently generated electric 
and magnetic fields. Plasma heating due to these shocks is 
the result of various microinstabilities that perturb the fields 
and particle distribution (Papadopoulos 1985, hereafter 
P85; Winske 1985). 

Collisionless shocks may be present when the Larmor 
radius is smaller than the central engine. The thickness of 
collisionless shocks is expected to be from several Larmor 
radii for quasi-perpendicular shocks to several hundred 
times the Larmor radius for quasi-parallel shocks (FKR92). 
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Note, however, that this result has not yet been confirmed 
for the near-relativistic plasma treated in this paper. For 
hydrogen-dominated plasmas the ion Larmor radius is 
roughly 

R Li ~ (0.03 cm )^~ l m~ lI2 M l 6 l2 r 314 



where @ = B/B v is the magnetic field in units of the virial 
magnetic field, which can be found from the relation 
Blj%n = GMpJR , and we have used equations (1) and (3). 
AGN central engines may support magnetic fields of up to 
the virial value of ~ 10 4 G at r ~ 20 (B90). The - 10" 8 G 
fields needed for the Larmor radius to be smaller than a 
perturbation in the central engine are easily obtained from 
matter accreted from stars. This limit is almost always satis- 
fied for the fields anticipated in the central engine (see, e.g., 
B90). 

For collisionless shocks in the central engines of AGNs 
one needs to use the kinetic equations for relativistic 
plasmas. Unfortunately, such equations are not fully devel- 
oped yet (Eilek & Hughes 1991), but several criteria from 
the studies of relativistic shocks can be used. For example, 
Barnes (1985) found that Alfvenic nonlinear waves steepen 
and become shocks. We therefore speculate that super- 
Alfvenic motions result in collisionless shocks. 

The Alfvenic velocity is 

C A = = ~ c$r~ 112 , ( 12 ) 

V 4™, m p 

which is comparable to c s for virial fields. The super- 
Alfvenic Mach number at that distance for Keplerian speed 
perturbations is thus 

= BJB . (13) 

It is therefore expected that strong collisionless shocks may 
form in the central engines for moderate magnetic fields, 
and moderate shocks may form for near virial 
fields. Instabilities that disrupt the thin disk may also 
result in the creation of a disk corona structure in which the 
corona is responsible for emitting the continuum for low 
enough accretion rates (Haardt & Maraschi 1991 ; Svensson 
& Zdziarski 1994; ZCC95). In the corona of a disk corona 
system the corona accretion rate is much smaller than in the 
tori models. In order to find the appropriate density for 
such corona one needs a model in which both disk and 
corona accretion are treated consistently. Such treatment 
was recently taken by ZCC95. Using Figure 5 in that paper 


for accretion rates of 

m ~ 0.001 

to m - 0.1 

and using 

equation (12) one finds 



X 

m 

Al 

< 

10 6 “ 7 cm s 

- 1 G~ l )B, 

(14) 

and thus at r ~ 20 





(3 x 10 3 “ 4 


(15) 


In both equations (13) and (15) the Alfvenic Mach 
number generally exceeds unity because the magnetic fields 
cannot quite reach their virial value. This is because energy 
from the magnetic field may be utilized for the outward 
transfer of angular momentum and/or for dissipation, 
which in collisionless plasmas may include, for example, 
reconnection events in the postshock region. 


3. FORMATION AND PROPERTIES OF THE SHOCKED 
MATTER 

3.1. The Shock Creating Perturbations 

Shock creating inhomogeneities in AGN central engines 
may be the result of nonlinear growth of global pertur- 
bation in thick disks that tend to form planetesimals 
(Papaloizou & Pringle 1984; Narayan 1990, and references 
therein), externally confined plasma sheets and filaments 
that cannot thermalize on orbital timescales (GR88; ST93; 
Celotti, Fabian, & Rees 1992, hereafter CFR92) or captured 
stars. The simplest physical shock creating perturbation 
considered in this paper is a main-sequence star. A more 
general treatment will be presented in our subsequent work 
(STC96; see also Caditz, Sivron, & Tsuruta 1996, hereafter 
CST96). It has been shown that an unbound star can be 
captured into a tightly bound orbit around a black hole 
through repeated interactions with a thin disk around the 
hole (SCR90), and that the probability of such an event in 
AGNs is even higher for stars in bound orbits (Pineault & 
Landry 1994). Once captured the orbit will be circularized 
and eventually the will be “ grounded down ” into the disk 
plane through repeated interactions with the disk. The time- 
scale for circularization can be shorter than that for plan- 
arization, which may happen very close to the central 
source. Therefore, in our model, the orbit of the captured 
star has been circularized, whereas the orbital plane is still 
inclined to the disk plane. This situation is expected for a 
wide range of the original relative inclination angles of these 
planes. 

The minimum conditions on perturbations with the 
potential to create shocks can be determined from the 
requirement that they accelerate to supersonic velocities. 
This condition is met when the gravitational force is greater 
than the drag force at , M = 1 : 

GMm ± , 

— (16) 

where p* ~ p g J( 2 cl is the pressure on the inhomogeneity 
surface (see § 3.2) and /), is the characteristic size of the 
inhomogeneity. This condition can be also written in the form 


P* 

f>9 


2 : 



(17) 


where p* is the average density of the inhomogeneity. This 
condition is easily met for a typical main-sequence star, 
which is used as an example shock producing inhomoge- 
neity below. 

Equation (17) may be even less restrictive if we include 
acceleration by radiation pressure and magnetic stresses. 
The amplitude and size of perturbations needed for shocks 
to form is therefore rather small, and shocks may form in 
the wake of large-amplitude nonaxisymmetrical waves in 
the disk. In this paper the possibility of nonlinear (soliton) 
shocks is not further developed because the dynamic evolu- 
tion of the perturbations is strongly dependent on dissi- 
pation (see STC96). 


3.2. Power Transferred to the Shock Front 

The power transferred from the star to the shock is at 
least 


^sts ~ P*)^* > 


(18) 
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where A* = 7 tZ> 2 is the cross section area of the star and p „ 
is the pressure near the stellar surface, behind the shock. 

The pressure at the leading edge of the stellar surface can 
be estimated from the treatment of shock formation in 
supersonic nonrelativistic flows past solids of revolution 
(Landau & Lifshitz 1987, hereafter LL87): 


P * > l-5pv\ 


[1 - 1/(5.// 2 )] 


1.5 


(19) 


The adiabatic constant was chosen to be its highest value of 
r = 5/3, which results in a lower limit on p*. The power 
transferred from the star to the shock is thus at least 


P sts - (1.5 x 10 43 ergs s l ) 


mMf, 


D*V 

Q / \^acc/ \ ) 


(20) 


The radiation thus emitted can reach the observed lumi- 
nosity with a single large perturbation, if the shocked 
material radiates efficiently. For collisionless shocks the 
ambient pressure is modified to be 


B 2 

Pa = p * + si’ 


(21) 


and the pressure from equation (19) is somewhat larger, 
especially with T ~ 4/3 for the heat capacity ratio when the 
magnetic effect is dominant. 

We next argue that if the accretion is already near its 
Eddington limit shocks will provide the mechanism for 
making the postshock region radiative. We separately treat 
the collisional and collisionless shock cases. 

3.3. Radiative Efficiency of Shocked Material 
3.3.1. Collisional Shocks 

In the near-Eddington regime radiative processes in the 
postshock region may be extremely efficient, owing to, for 
example, bremsstrahlung from thermal pairs. We speculate 
that most of the energy is radiated by pairs that are created 
in the postshock region. The presence of some pairs is 
expected in collisional hydrogen plasmas with bulk motion 
with velocity of order 



(see 189) and in the central engine the bulk motion is of order 
v+, which is generally much larger than U in equation (22). 

A more significant pair creation is present in the post- 
shock region because the temperature of the postshock elec- 
trons is & e > 1, and because these electrons are 
thermalized. We hereafter denote with a prime the post- 
shock quantities. In this scenario the optical depth for pair 
annihilation is larger than unity, and a high-density pair 
equilibrium is established in the postshock region. These 
circumstances were investigated by Done & Fabian (1989, 
hereafter DF89), CB90, and others. In their simulations 
these authors found that the spectrum, which is a power law 
due to, for example, Comptonization of bremsstrahlung 
radiation, is flattened because of the increased compactness 
of the source. 

This situation is also likely to apply locally to the post- 
shock region because the optical depth for pair annihilation 
is of order V - 0An\ <r T A = 5© 2 and © e > 1. The timescale 
for mildly relativistic (0, - 1) electrons to thermalize t e ~ 
l/(ca T n-ln A) is shorter than the shock crossing timescale 
f sw ~ X/c s , and, using equation (2), we find the following 


ratio to be smaller than unity: 



3 c s 1 

— < 1 

64 c © 2 ~ 


Here we have used equation (5) and the condition 


(23) 


^ :> 4 , (24) 

from 189. The condition of equation (23) is generally ful- 
filled, and the postshock electrons are therefore thermalized 
(see also CST96 and STC96, for a more detailed approach). 

For the thermalized postshock pair equilibrium 
described above, the total number density of pairs, n\, may 
be estimated by evaluating the postshock energy of the ions 
and assuming equipartition. The ion energy is approx- 
imately 



where we have used the postshock energy from 189 and 
equation (24). 

At r ~ 20 the bulk velocity is roughly u* ~ c/20 1 ' , and 
the postshock temperature is roughly 0' = 30 mjm p , which 
translates to more than ~40 electrons and positrons per 
postshock ion in equation (25), and we have used equations 
(8) and (22). Our assumption that at the above location the 
postshock creation of pairs is unstable for m ~ 0. 1 is sup- 
ported by Bjomsson & Svennson 1991 (see also Fig. 1). 

Because the density of pairs in the postshock region is so 
much larger than that of ions in the preshock region, and 
because bremsstrahlung radiative emmissivity is pro- 
portional to the square of the density, the postshock pairs 
may lose most of their energy to bremsstrahlung radiation. 
The total bremsstrahlung power per unit volume is of order 

j > (1.2 x 10 -23 ) erg cm -3 s _1 (n'ff(&' e ) 112 , (26) 

(e.g., FKR92). The volume of emitting postshock matter is 
of order X tiD 2 , and the power radiated assuming an opti- 
cally thin source is 

P b ~ jXnDl , (27) 

which, using equations (20), (25), and (26), generally exceeds 
P sts . The power in equation (20) is therefore emitted as 
bremsstrahlung radiation if P b > P sts , which is equivalent to 
the condition 

(0' e ) ,/2 0 2 ^O.lr-‘ , (28) 

which is usually satisfied, and we have used equations (2), 
(3), and (20). In Figure 1 we show that the increase in tem- 
perature and electron-positron number density across a 
typical central-engine collisional shock results in an 
increase by several orders of magnitude in bremsstrahlung 
emission rate across shocks of moderate to high Mach 
numbers. The Rankine-Hugoniot relations are taken from 
189, the density of pairs in both the preshock and postshock 
regions is taken from Svensson (1982, 1984) and Bjornsson 
& Svensson (1991), and the typical parameters of the two- 
temperature plasma are taken from WL89 at a distance of 
r ~ 20. This figure is taken from a more general discussion 
in CST96. Figure 1 gives the lower limit on the postshock 
radiation, because additional increase in emission due to, 
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M 

Fig. 1. — The effects of the increase in temperature, number density, and number of pairs across a typical central engine collisional shock with Mach 
number M results in an increase in bremsstrahlung emission rate. The postshock-to-preshock ion temperature, ion density, and bremsstrahlung emissivity 
ratios, T'/T , N'/N, and b’/b, are the dashed, dash-dotted, and solid curves, respectively. They are plotted vs. the Mach number. To the right of the thick solid 
line is the pair runaway region where condition (25) is fully satisfied. The pair rate equations the preshock and postshock regions in our code i are taken from 
Svensson (1982, 1984) and Bjornsson & Svensson (1991). The Rankine-Hugoniot relations for relativistic collisional plasmas are taken from 189. This figure is 
taken from a more general discussion in CST96. The preshock temperatures are taken from WL89 for m ~ 0.1 at r ~ 20. The postshock ion temperature/ 
electron temperature ratio of roughly ~ 10, is appropriate for log N'/N ~ 0.6. 


for example, Comptonization of the bremsstrahlung radi- 
ation, for example, is not taken into account. We therefore 
conclude that the postshock region can easily radiate away 
most of the power transferred from the star to the shock. 
We use this result in § 4 for deriving a possible light curve 
from an emitting shock. 

It has been suggested that for certain ranges of com- 
pactness the hot gas in the central engine is unstable, and 
the perturbations may result in the formation of a “ two- 
phase” region consisting of a dense cool (not completely 
ionized) phase and a hot phase (GR88). In a similar manner 
the postshock matter in our model may also form a two- 
phase configuration. The postshock compactness is of order 



where we have used equations (3) and (20) and assumed 
<5 — 0.01. This compactness is related to the standard com- 
pactness, / = L<r T /(Rm e c 3 ) (Guilbert, Fabian, & Rees 1983). 
As in GR88, the high compactness leads to an interesting 
effect: In some postshock regions ions are trapped and may 
lose most of their energy to the surrounding electron- 
positron pairs through a runaway cooling, which involves 
Comptonized bremsstrahlung radiation and then bound- 
free transitions. Counter to our experience with other 
shocks some postshock ions after being initially heated may 
cool below the original preshock temperature to a typical 


GR88 cold phase temperature of 

T L l 1 ' 4 

r - - L (1 " A) > (30) 

where T ps is the postshock temperature and A is the flux 
averaged albedo. This matter, which probably forms clouds 
of denser cold matter, may be “warm,” i.e., partially 
ionized, because of its close proximity to the strong radi- 
ation from the nearby shock, or “ cold ” if the lifetime of the 
cold matter in the cloud is longer than the Keplerian time- 
scale and the matter is not intensely irradiated. 

A configuration with cold plasma in the central engine 
may therefore exist. The cold component may give rise to 
some emission and absorption features such as the gravita- 
tionally redshifted Fe K line recently observed by ASCA 
(Tanaka et al. 1995). The cold (or warm) phase may also 
eclipse the hot plasma under certain conditions. The eclipse 
may modify the primary emission by absorption in a post- 
shock foamy sheet behind the Mach surface if P sts is 
expended on the cooling of the postshock matter, and if the 
lifetime of the cold blobs is long enough (see Appendix). In 
subsequent papers (CST96; STC96) we will present a full 
numerical solution for the rate equations in the two-phase 
postshock region, which accounts for these effects. 

3.3.2. Collisionless Shocks 

A naive application of the same procedure used in § 3.3.1 
to collisionless shocks would have us compare the shock 
crossing time, which is either t sw ~ r Li /c A or the possibly 
longer t sw ~ (r Li /c A )( 1 + a> b r Li /c A ), which is increased due to 
confinement (where co b = eB/mc is the cyclotron frequency), 
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to the deflection time X/(2KT/m) 112 (which is the typical 
timescale for, e.g., conduction losses). Shkarofsky, Johnston, 
& Bachynski (1966, p. 23) used these timescales t sw for the 
diffusion time out from a confined region. The problem with 
this approach is twofold: (1) The width of a collisionless 
shock tends to be at least several Larmor radii (e.g., 
KEH85) and the crossing time is therefore several times t B . 
(2) As discussed in § 2.2 the plasma losses across collision- 
less shocks are determined not by single binary collisions, 
but rather by collective interaction between particles and 
self-consistently generated electric and magnetic fields. 
Most of the losses in collisionless shocks, (especially in the 
near virial fields limit), may be due to, for example, various 
microinstabilities in the collisionless plasma rather than 
conduction (Winske 1985; P85). 

At present there are no analytical Rankine-Hugoniot 
relations for relativistic collisionless shocks in matter with 
density in excess of n ~ 10 10 cm -3 . We therefore param- 
eterize the increase in magnetic field magnitude, and poss- 
ibly particle density and velocity, across the shock. This 
increase would result in an increase of the synchrotron and 
synchrotron-self-Compton radiation and would take place 
over a timescale similar to t B . 

We speculate that the Rankine-Hugoniot relations for 
the conditions envisioned in the central engine have some 
common features with the Rankine-Hugoniot relations of 
lower energy collisionless shocks (e.g., P85, and references 
therein). For simplicity’s sake we assume that the quiescent 
radiation is due to synchrotron radiation only, such that 

L(0)~2 x 10 ~ 9 B 2 ynV x , (31) 

as in B90, where L(0) is the actual observed bolometric 
luminosity and V x is the volume of nonthermal emitting 
medium, which can be either in a quasi-spherical form or in 
a corona-over-a-disk form (see § 2). The resulting postshock 
radiation density is therefore of order 





By uo) 

B 


(32) 


We hence estimate the overall compactness of the post- 
shock matter for quasi-spherical accretion : 


strong shock with high Jl A number (see § 2.2) the postshock 
parameters may result in a significant increase of the flux 
and the steepening of the spectrum. We discuss the implica- 
tion of this result in STC96. 

4. MODIFICATION OF THE CONTINUUM SOURCE 
4. 1 . Effects of Radiative Shocks 

Since the speed of the “star” in our model is assumed to 
be roughly ~0.2c at r = 20 we will neglect the increase in 
shock radiation due to the Doppler effect (for corrections 
see STC96). We note, however, that the angle spanned by 
the shock front area vector and the line of sight determines 
the percentage of outcoming shock radiation due to the 
postshock eclipsing, which is dealt with in the next sub- 
section. Therefore, emission only doubles when the shock 
front is directed toward the observer (see Fig. 2). In our 
approximation we assume that the spectral index in the 
postshock pair-dominated region is ~ 1, as was shown by 
DF89, CB91, TK95, and others, whereas in the ambient 
matter the spectral index is —0.5, as is expected in synchro- 
tron and synchrotron-self-Compton sources (B90, and ref- 
erences therein). 

In Figure 3 we show the effect of a shock on the spectral 
index in the soft X-ray waveband in the 0.2-2 keV window, 
the hard X-ray waveband in the 10.0-20.0 keV window, and 
the very hard 100.0-300.0 keV window. The angle /? is 
spanned by the shocking star, the central engine and the 
observer-black hole line of sight. One can clearly see that 
the outcoming flare is soft as observed in several compact 
Seyfert galaxies. The lesser softening in the higher energy 
bands is also consistent with observations (see MDP93). We 
assumed that the maximum observed flux with spectral 
index of order unity in the 0.1-300 keV range is similar to 
the flux from the ambient matter. We note that a spectral 
break near 300 keV was assumed in our work because we 
assumed a two-dimensional slab corona configuration for 
our pair spectrum, as in TK95, which results in a spectral 
break near ~300 keV for high enough compactness. A 
combination of observations with present and future mis- 
sions, such as XTE , IMM, and AXAF, may be able to 
determine if the light curves fit these characteristics. 


Ps^Ri 


= 10 m: 


m e c" 

and for coronal accretion : 


.(JL.V.1) 

\0.0i;\20/ 



, (33) 


P 3 <t t Rx h 

m e c 3 



(34) 


Extrapolating from interplanetary collisionless shocks 
(KEH85) and extragalactic collisionless shocks in jets (see 
Eilek & Hughes 1991, and references therein) we speculate 
that the postshock magnetic field B’ = kB and density n' are 
larger than the corresponding preshock parameters. The 
compactness may therefore be increased enough for pair 
production to be important. The resultant pair production 
will be nonthermal, and the steepening of the spectrum of 
emission from the postshock region may be significant. This 
postshock emission will result in steepening of the overall 
spectrum if the postshock emission is comparable to that of 
the quiescent emission. Given a large perturbation and a 


4.2. Effects of Postshock Cold Matter 

Here we obtain the light curve for radiation from hot 
plasmas in the central engine, which is “eclipsed” by the 
cold postshock matter. The central engine will be eclipsed 
when the postshock cold matter creation timescale, of order 
t e in equation (23), is shorter than the cold matter timescale 
t s (eq. [A9] in the Appendix), and when the postshock 
matter column density is of order n'c s t s ~ 10 21-24 cm -2 
(depending on the ionization state of the postshock “ cold ” 
matter). Both conditions are met for postshock matter in at 
least some AGNs (see eqs. [2], [3], [7], [24], and [A9]), and 
an eclipse of the central engine is therefore possible. 

In obtaining the light curve our approach is geometrical 
in nature. We assume that a conical shock front, similar in 
shape to a simple Mach cone (see § 2) eclipses a circular 
continuum source. To highlight the geometrical effects we 
also chose a rectangular continuum source. A possible 
picture of what happens near the central region of an AGN 
with an inhomogeneity is shown in Figure 2. Here a shock 
inducing star (at the tip of the cone) orbits inside a thick 
accretion torus, which is located outside the X-ray contin- 
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uum region but confined to a region smaller than the thin 
outer disk (points FI and F2 are the inner edges of the thin 
disk). From an observer’s point of view the points O, P, and 
Q on the Mach cone (Fig. 4) form a “ fan ” in the shape of a 
triangle that partially obscures the internal continuum 
source. The type of occultation depends on the optical 
depth of the matter on and near the shocked surface. 

Since the extent of the optically thick shock may be 
limited, as was shown in § 3, we choose the point of view of 



Fig. 3.— Effect of a shock on the spectral index of flux from the central 
engine with a radiative shock front is plotted vs. the angle spanned by the 
Keplerian perturbation, the central black hole and the line of sight The 
solid line is for flux in the “ soft " 0.2-2 keV window, the dashed line is for 
“ hard ” 10 ~ 20 keV window, and the dash-dotted line is for the 
100-300 keV window. The trend of harder flares in the higher energy 
bands is apparent. At an angle n/2 the Mach cone is pointing toward the 
observer. 


an observer in direction G, which requires the shortest 
extent of optically thick matter. Figure 4 is a simplification 
of the observer’s view, which is used for the calculation of 
the light curve. It is based on the following assumptions : 

1. The continuum X-ray emitting region is assumed to be 
circular, with homogeneous surface brightness, as viewed 
when unocculted. Its radius is R = We note that in an 
alternative model, in which the soft X-ray region is larger 
than the hard X-ray region, the soft minimum flux will be 
significantly delayed with respect to the hard flux. 

2. The Mach surface is assumed to be an equilateral tri- 
angle, passing across the center of the circular source. 

3. E>ensity variations are simplified by assuming three 
density regions. In region 1 (angles smaller than 0,) the 
matter is opaque to both hard and soft X-rays. In region 2, 
with angles between 0, and 0, the matter is transparent to 
hard X-rays, but opaque to soft X-rays. In region 3, behind 
the darker triangle, the density of matter varies with dis- 
tance, due to the spread of matter, and is partially transpar- 
ent. We consider only the hard X-ray waveband in the 
10—20 keV window in which no absorption is expected and 
the soft X-ray waveband in the 0.5-2 keV window in which 




Fig. 4.— Observer’s view of the AGN in our model with covering factor 
off = 0-5. This is a simplified form of the view from point G in Fig 2 The 
stepwise fashion, which was chosen to describe the absorber, is for demon- 
strative purposes. See § 4 for details. 
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F,g. 5.— Light curves near a dip that will be observed when the source 
is eclipsed by the postshock matter is shown for the softer (dash-dotted 
curve) and the harder ( solid curve) X-ray windows with <j> - 45. The short 
dashed curve is for the case of a smaller Mach surface opening angle of 
<b = 30° for the hard X-rays. The long-dashed curve is for the simplified 
case in which a triangle with an opening angle <j> = 45° eclipses a square 
source. In all cases the maximum covering factor is/ - 0.5 (see § 4 tor 
details). 


absorption due to cold or warm oxygen is expected^ We 
chose the maximum covering factor to be /- 0.5 lor 
demonstrative purposes. 

4. The strength of the predip flare was assumed to be 
negligible for demonstrative purposes. 

In Figure 5 we see the light curves near a dip in two 
energy bands that will be observed when the source is 
eclipsed by the postshock matter. With a Mach surface 
opening angle of <f> = 45°, the softer X-ray light curve is 
represented by the dash-dotted line and the harder X-ray 
light curve by the solid line. The dotted curve is for the case 
of a smaller opening angle of (j> = 30°, for the hard X-rays. 
The dashed curve is for the simplified case m which a tri- 
angle with an opening angle 4> = 45° eclipses a square 
source. In all cases the maximum covering factor is /= 0.5. 

The decrease of flux toward the dip is energy dependent. 
A delay is observed in the sense that the decrease in photon 
count in higher energies lags the decrease in the lower ener- 
gies. The light curve in the hardest window has a cusphke 
shape, dropping to half its value at the bottom. The bottom 
was chosen to be at h alf the intensity so that the size of the 
central sector of absorbing material within </>!, which is 
opaque to the harder X-rays, is approximately half the size 
of the harder X-ray source size during maximal occultation. 
The maximum value for R x is thus obtained with <P i - 45 • 


5. DISCUSSION 


The results of §§ 2, 3, 
bation, an accretion rate 


and 4 imply that given a pertur- 
and a black hole mass one could 


calculate temporal features and confirm the 
“ condensation ” of cold (not completely ionized) matter. 
Our work thus complements the work by GR88 and LW88. 
We also argue that one could also calculate the variability 
and power spectrum, the rate of angular momentum trans- 
fer, and the effects of particle acceleration. The results of §§ 2 
and 3 may give rise to the prediction of the r hs = ljl s ratio, l h 
and l s , where l h is the compactness of the hard nonthermal 
radiation and l s is the compactness of the soft thermal “ seed 
photons ” (see STC96). 

When the accretion rate is relatively low one expects a 
disk corona configuration. In this configuration accretion 
disk instabilities are expected to yield microvariability (e.g., 
\Y93). argued in §§ 2 and 3 collisionless shocks will then 
tend to spread these perturbations throughout the corona 
at the fastest magnetosonic speed, thus enabling large 
amplitude variability. Since the coronal compactness is 
relatively small only a fraction of the perturbation energy 
will be radiated away (and most of the perturbation may 
travel outward as winds). Note, however, that various 
microinstabilities may increase the efficiency of energy 

transfer. . . , , . • 

The radiation will tend to be relatively hard since no 
Comptonized bremsstrahlung radiation is expected and 
only synchrotron-self-Compton (SSC) or nonthermal pairs 
cascades are possible locally, with a very hard outcoming 
spectrum (e.g., B90; TK95). The compactness l s also 
depends on m (FKR92), and r hs may not change as fast as (*. 
The radiation may be SSC only if l h + l, is such that small 
IR blocking cold cloudlets form in the shocked material (as 
in CFR92) or else TR- X-rays correlation should be observ- 
able (Done et al. 1991). Reflection of the hard radiation 
from the relatively cold disk component may result in an Fe 
emission line, such as the one recently observed by ASCA 
(Tanaka et al. 1995, and references therein). 

When the accretion rate is relatively high (m > 1) one 
may expect a very dense corona or a thick two-temperature 
torus. Perturbations will be the result of instabilities in the 
disk (e.g., Narayan 1990, and references therein) and the 
resultant shocks will produce bremsstrahlung-Compton 
radiation or may even result in a two-component Guilbert 
& Rees medium. The emission from such media becomes 
softer when l h + l s become higher (ST93; Nandra & George 
1994). In general one expects l s to increase more than l h in 
this case, because of the soft photons that originate or are 
reprocessed in the cold component. 

The emission spectrum for the high m is dominated by a 
more thermal-like component with a possible high- 
temperature bremsstrahlung, which may imply AGN evolu- 
tion, which is thought to be responsible for the observed 
X-ray background. 

Shocks can naturally explain a common feature of the 
spectra of most AGNs. In the spectra of radio-quiet AGNs 
the ratio of thermal to nonthermal emission is probably of 
the order 0.3 < L T jL NT < 3, where h T is the thermal lumi- 
nosity and L nt is the power-law nonthermal luminosity 
(LW88). This result is readily explained if the hot contin- 
uum emitting plasma is nearly optically thick (Sivron 1995). 
A mechanism in which the optical depth is limited to unity 
involves a two-phase medium : A medium that would have 
been optically thick if the accretion were spherical becomes 
a two-phase medium with decreased optical depth for the 
hot phase whenever the optical thickness exceeds unity, or 
equivalently whenever the luminosity is near-Eddington 
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(ST93; STC96; GR88). Such two-phase media are naturally 
produced with shocks, because shocks are so readily avail- 
able in the central engine (as was shown in § 2) and because 
shocks naturally channel the energy of perturbations into 
radiative cooling, which cools portions of the hot phase in 
the case of near-Eddington luminosities (as was shown in § 
3). The cycle is completed when the speed of sound in the 
rarefied hot phase again falls below the Keplerian velocity 
due to the “ evaporation ” of cold phase clouds, which is 
usually faster than the dynamical timescales, or due to 
cooling processes in the hot plasma, and shocks are again 
possible (for a typical timescale of “ evaporation ” see, e.g., 
Zel’dovich & Raizer 1967, p. 571 and the Appendix). 

If radio-quiet AGN central engines have a disk corona 
geometry, the observed spectra cannot be explained by 
homogeneous coronas (Fabian et al. 1995; Stern et al. 1995; 
Iwasawa et al. 1996). This paper establishes the presence of 
shocks that naturally disrupt the homogeneous corona. 
Perturbations in the hot homogeneous coronas may be 
considered in the sense described in this paper in order to 
establish the true geometrical shape of the inhomogeneous 
corona and calculate the resultant spectrum. 

Shocks may also explain the observed power spectrum of 
the central engine. A weakness of our assumption that a 
single star is responsible for the shock and resultant emis- 
sion is that periodicity, which has not been observed in 
AGNs, may be expected. A star orbiting a supermassive 
black hole can be a curiosity, or a rule — depending on 
future observations, but from § 3.2 it is apparent that any 
large perturbation present in the innermost 20 R g can yield 
the luminosity change of equation (20). It is therefore prob- 
ably more likely that the power spectrum is the result of 
perturbations that are not associated with a captured star. 
If the conditions in equations (4) and (30) are met for these 
perturbations, and if perturbations of characteristic size, 
time, and brightness occur at a random distance from the 
black hole, the resultant power spectrum slope will be 
between that of random shot noise and chaotic noise. Our 
subsequent work (STC96) uses the arguments of Sunyaev & 
Titarchuk (1980) and DF89 to obtain these power spectra, 
which are comparable to the power spectra observed in 
radio-quiet AGNs (Green 1993, and references therein). The 
condition in equation (30) on the compactness then deter- 
mines that a one-dimensional introduction of pertur- 
bations, rather than the two-dimensional introduction of 
perturbations as in thin disk models, is needed in order to 
reproduce the observed power spectra. The physics of 
random introduction of perturbations in one dimension 
may be easier to explain than the random introduction of 
perturbations in two dimensions (STC96).Condition (29), or 
an equivalent for a nonspherical perturbation, may result in 
a high-frequency cutoff in the power spectrum (STC96). 

If the shocks are collisionless the source may not be opti- 
cally thick, and modification through pair creation may not 
be as significant, unless the fields are near virial. In the case 
of collisionless shocks that do not satisfy the conditions in 
equations (32)-(34), one should note that energy may still be 
deposited by the shock to the ambient hot plasma. Various 
microinstabilities may still result in microvariability, but 
large-amplitude variability is harder to produce. 

In the above discussion of variability we did not try to 
calculate the timescale in which the postshock matter 
returns to its original state. We generally assumed that the 
recovery timescale is shorter than the dynamical timescales. 


If the recovery timescale is larger than the dynamical time- 
scale, the power spectrum is further modified. The effect of a 
two-phase medium on the variability is also calculated in 
our subsequent work (STC96). 

The spectral variability of flares and dip described in §§ 3 
and 4 was already found to be similar to that of several 
Seyfert 1 galaxies such as NGC 4051 (Kunieda et al. 1990), 
MCG — 6 — 30 — 15 (Nandra, Pounds, & Stewart 1990; 
Matsuoka et al. 1990; Fabian et al. 1994; Otani et al. 1995; 
Sivron et al. 1996a), and others (MDP93). The similarity is 
not only in the sense that the power-law spectral slope is 
anticorrelated with the overall flux, but also in the sense 
that the low flux level is correlated with some emission and 
absorption features. Although this correlation was tradi- 
tionally explained to be the result of matter in the line of 
sight, which is far from the central engine (e.g., Nandra et al. 
1990) it is generally hard to explain why this matter is 
always smaller in thickness than the size of the central 
engine (Sivron et al. 1996a). The emission and absorption 
features may thus be associated with an increase in optical 
depth of cold matter within, or very close to, the central 
engine, as in the case described in § 4.2. 

Shocks and their by-products may be used for the 
outward transfer of angular momentum in the central 
engine (STC96). Spiral shocks, which were considered for 
the outward transfer of angular momentum in the outer 
disk, may form when a massive companion perturbs the 
disk (Savonije, Papaloizou, & Lin 1994, and references 
therein, in reference to galactic sources). Angular momen- 
tum can be transferred in this way in AGN inner disks as 
well (Chakrabarti 1995). Spiral shocks probably take a 
completely different geometrical shape in compact central 
engines, where the geometry of the emitting material is not 
that of a thin disk. The angular momentum transfer in the 
central engine in high-accretion Seyfert 1 galaxies and 
radio-quiet quasars can be due to chaotic, rather than 
spiral, shock structure (STC96). 

In the case of lower accretion rates significant outflows 
must also be considered as a part of the total angular 
momentum transfer, since the optical depth is limited and 
mass ejection may become significant. However, the avail- 
ability of shocks in the central engine and the presence of 
mechanisms that can disturb the axisymmetrical structures 
makes it imperative that shock-related angular momentum 
transfer is significant. Examples of mechanisms that can 
disturb the axisymmetry include a poloidal magnetic field 
structure, a Lense-Thirring precession of tilted axisym- 
metrical structures in a Kerr geometry, etc. 

Shocks may also be an indirect agent in the outward 
transfer of angular momentum. The postshock cold clouds 
may transfer angular momentum through binary and triple 
collisions (STC96). Because the hot component in a GR88 
two-phase model is fairly optically thin, some cold clouds 
may also be ejected carrying with them large amounts of 
angular momentum. 

Other processes affected by the presence of shocks 
include particle acceleration and pair production. Many 
models use an input energy spectrum for a large nonthermal 
electron population in order to obtain their model spectra 
(B90, and references therein; CB90). From our shock model 
one should directly calculate the normalization of this dis- 
tribution by utilizing, for example, the theory of Fermi and 
shock-Fermi acceleration for a specific postshock preshock 
compression ratio (Blandford & Eichler 1987; Levinson 
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1994). We note that modifications due to Fermi acceleration 
will probably result in the lowering of cooling efficiency, as 
in equation (30). That a significant portion of the particles 
may gain a nonthermal high-energy distribution function 
was anticipated (but not calculated) by many previous 
authors that discussed pair modifications of spectrum and 
structure of the central engine (see Guilbert et al. 1983; 
Begelman, Sikora, & Rees 1987; Tritz & Tsuruta 1989; 
Kusunose & Takahara 1988; CB90; GH94). 

6. CONCLUSIONS 

We have shown that shocks in the central engines of 
AGNs form for a wide range of accretion parameters. We 


have also shown that perturbations in the central engine 
will effectively radiate their energy away through shocks 
and that for a wide range of accretion parameters that radi- 
ation will be observable. Given the generality of our results 
we believe that shocks should be fairly common and that 
their effects on observations should be significant. 

We thank the anonymous referee for some very helpful 
remarks. We thank M. J. Rees, N. Iwamoto, and N. Sivron 
for useful discussions and comments. We also thank W. 
McHargue and M. J. Kellen for their help with supporting 
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APPENDIX 


The lifetime of the thermal gas sheets and filaments that form behind the shock front can be evaluated using one of the two 
following methods : 

The minimum lifetime of plasma sheets and filaments can easily be approximated for the case in which the sheets are in 
the shape of a spherical cloud and all of the incident radiation is transformed into an energy increase of size A E : 


t, > ^ = 140r 1/2 
F 


\cJ\tim H c 2 J ’ 


(Al) 


where F is the maximum flux incident on a cloud, r c is the radius of a cloud and T c ~ 10 8 K is the final temperature of the 
clouds in which all the elements are completely ionized (see also Bond & Matsuoka 1993). For a typical cloud of size 10 8 cm 
(for M ~ 10M 6 , using the same arguments as in CFR92), the lifetime of the clouds is greater than 10 4 s. This result is a lower 
limit, because no radiative losses are considered in spite of the fact that bound-free and free-free emissions, which are the 
dominant mechanisms in the creation of such thermal matter, have already been proven to be efficiently radiating energy in 
§ 3.2. Whereas this timescale describes the dissipation of a cloud, we next try to determine the beginning of such a process, 
which we show to have a longer timescale. 

We assume that the temperature decrease is due to adiabatic expansion into vacuum and calculate the subsequent timescale 
necessary for the optical depth to be halved. The following equations govern the adiabatic expansion of clouds. The 
continuity equation in spherical coordinates, for a spherically symmetric time-dependent case (LL87), is 



\ T ( r V) = 0 > 


(A2) 


where r is the radius of the exploding plasma cloud and v is the velocity of a fluid element. The Euler equation is 


dv 1 dP 
dt + p dr 


(A3) 


The adiabatic (energy) equation is 


p = 7’y/(y-i) > (A4) 

where y = 5/3 is roughly the adiabatic index for an ideal monoatomic gas, for hydrogen-dominated plasmas. The equation of 
state is 


The total mass of a cloud, M cld , is 



47t , 

Mm = yP' 


(A5) 

(A6) 


The optical depth of a cloud, assumed here to be bremsstrahlung and bound-free dominated, is 

t ~ 1.0 X 10 3 rT~ U2 Z 2 p 2 (~^ (A7) 

The expansion of the clouds results in their transparency. The optical depth of the clouds as they expand (and are assumed 
to be opaque, t = 1, at that size), can be related to the corresponding variations in cloud size, density, and temperature 
without the need to solve equations (A2) and (A3) explicitly. By substituting equation (A6) into equation (A5) and then into 
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equation (A4), we get T oc r 3 3y , and substituting the result into equation (A7) we find 


*(r 0 ) 



(A8) 


Here r f is the final size of the expanding cloud. From this we can see that a bremsstrahlung and bound-free dominated, 
adiabatically expanding cloud will only need to grow by r f /r 0 ~ 1.12 for the optical depth to decrease by a factor of 2. 
Substituting these initial and final sizes into equation (A3) we can estimate the lifetime of a cloud. Since the pressure, density, 
and size change very little when a cloud opacity is halved, the Euler equation becomes roughly d 2 rjdt 2 ~ — (VPoX^o/ r o)- 
Integrating twice, using the initial conditions v(t = 0) = 0 and r(t = 0) = r 0 , we get 


t s > -(100s) 



(A9) 


where r 8 = r o /10 8 cm, T 6 = T/10 6 K, r 0 andP 0 are the initial size and pressure of the dense sheets, and is the time from the 
beginning of expansion. Hence, the postshock sheets become transparent on a timescale of hundreds of seconds. As a result a 
“curtain” of cold matter of thickness t s c s is comoving with the Mach cone. This curtain may partially obscure the radiation 
from the central engine if its optical depth, /i7 s c s cr, where ri is the postshock cold matter average density and a is the cold 
matter’s effective cross section, is nearly unity. Because of confining pressure from magnetic fields, which may appear on the 
right-hand side of equation (A3), the lifetime of a cloud may be longer than that in equation (A9). 

Considering the velocity of the inhomogeneity found in § 2 and the fact that the postshock matter moves at subsonic 
velocity, this result determines that the clouds form a “foamy” wall of thickness ~c s t s ~ 10 10 cm. A typical filling factor of 
C ~ 0.01 is reasonable for the ratio of ambient and thermal matter under the assumption that half the ambient matter 
becomes thermal in the postshock region (with the necessary covering factor ~ 1 ; see CFR92 and ST93). Therefore, on the 
average, one sheet of thickness 10 8 cm can partially eclipse the continuum source. 
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